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ABSTRACT 

Context. Age and mass determinations for isolated stellar objects remain model-dependent. While stellar interior and atmospheric 
theoretical models are rapidly evolving, we need a powerful tool to test them. Open clusters are good candidates for this role. 

Aims. We aim to create a fiducial sequence of stellar objects for testing stellar and atmospheric models. 

Methods. We complement previous studies on the Hyades multiplicity hy Lucky Imaging observations with the AstraLux Norte 
camera. This allows us to exclude possible binary and multiple systems with companions outside 2-7 AU separation and to create 
a "single-star sequence" for the Hyades. The sequence encompasses 250 main-sequence stars ranging from A5V to M6V. Using the 
"Tool for Astrophysical Data Analysis" (TA-DA), we create various theoretical isochrones applying different combinations of interior 
and atmospheric models. We compare the isochrones with the observed Hyades single-star sequence on J vs. J - Ks, J vs. J - H and 
Ks vs. H — Ks color-magnitude diagrams. As a reference we also compute absolute fluxes and magnitudes for all stars from X-ray to 
mid-infrared based on photometric measurements available in the literaturefROSAT X-ray, GALEX UV, APASS gri, 2MASS JHK,, 
and WISE W1 to IV4). 

Results. We find that combinations of both PISA and DARTMOUTH stellar interior models with BT-Settl 2010 atmospheric models 
describe the observed sequence well. We use PISA in combination with BT-Settl 2010 models to derive theoretical predictions for 
physical parameters (Teff, mass, log g) of 250 single stars in the Hyades. The full sequence covers the mass range 0.13 to 2.3 Msun, 
and effective temperatures between 3060 K and 8200 K. 

Conclusions. Within the measurement uncertainties, the current generation of models agree well with the single-star sequence. The 
primary limitations are the uncertainties in the measurement of the distance to individual Hyades members, and uncertainties in the 
photometry. Gaia parallaxes, photometry and spectroscopy will greatly reduce the uncertainties in particular at the lowest mass range, 
and will enable us to test model predictions with greater confidence. Additionally, a small (~0.05 mag) systematic offset can be noted 
on 7 vs. 7 - X and K vs. H - K diagrams - the observed sequence is shifted to redder colors from the theoretical predictions. 

Key words. Techniques: high angular resolution - Stars: atmospheres - Stars: binaries - Stars: distances - Stars: fundamental 
parameters - Stars: general - (Galaxy:) open clusters and associations: individual: the Hyades 


1. Introduction 

The knowledge of physical parameters of stellar and sub-stellar 
objects provides the basics for astrophysics. Observations allow 
us to determine distance (using parallaxes), age (through mem¬ 
bership in clusters or moving groups), mass (by calculating or¬ 
bital parameters in binary systems), and radius (in transiting and 
eclipsing systems). Asteroseismic scaling relations can also pro¬ 
vide estimates for stellar mass and radius. However, these rela¬ 
tions mays suffer from unknown systematics effect and require 

* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Institute de Astroflsica de An- 
dalucfa (CSIC) 


additional calibration using observations and theoretical models 
(e.g. |Miglio et al.|2012| l. Therefore, when dealing with isolated 
stellar objects, we usually have to rely on theoretical models to 
determine mass and age. While stellar and atmospheric models 
are rapidly evolving, we need a powerful tool to test and calibrate 
them. Open clusters are good candidates for this role, since they 
contain many coeval objects of the same chemical composition 
spanning a range of masses, thus avoiding the problems of small 
number statistics. Additionally, open clusters have distance and 
age estimates that are independent of theoretical models to be 
tested. 

In general, the validation of theoretical models consists of 
two independent problems - one is to test stellar interior struc- 
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tures, another is to inspect atmospheric models. However, this 
process is non-trivial because interior structure models provide 
us with physical parameters, such as effective temperature, lu¬ 
minosity, mass and age, that need to be converted into observed 
quantities (magnitudes, colors and fluxes). This conversion ex¬ 
ploits grids of synthetic spectra provided by atmospheric models 
and requires careful calibration for specific photometric bands 
(e.g. |Da Rio & Robberto|2012| l. 

In the past, open clusters have already been used as bench¬ 
marks for theoretical models. [Bell et al.| ( |2012| l used the Pleiades 
to test existing pre-main sequence isochrones by comparing their 
predictions with well-calibrated color-magnitude diagrams in 
the wavelength range of 0.4-2.5 pm. Bell et al. have shown that 
no pre-main sequence model can describe the observed Pleiades 
sequence for the temperatures cooler than 4000 K. The predicted 
fluxes are over-estimated by a factor of 2 at 0.5 pm, with the dif¬ 
ference decreasing with increasing wavelength. 

In comparison to the Pleiades (d ~ 120-140 pc; Percival et al. 
[2005 1 [van Leeuwen 2009|l, the Hyades open cluster is closer to 


the Sun (d ~ 45 pc; Perryman et al.|19^8j van Leeuwen]200^ . 

This allows us to resolve companions with smaller physical sep¬ 
arations, and also to analyze less luminous, hence less massive, 
objects. Moreover, [Roser et al. ( |201 Ij ) reported 724 likely mem¬ 
bers of the Hyades, with individual kinematic distance estimates 
using the convergent point method (e.g. [van Leeuwen|2009| l. In¬ 
dividual distance measurements allow us to get more precise ab¬ 
solute magnitudes for each member of the Hyades. The previous 
attempt to test stellar models with the Hyades ( [Castellani et al.j 
[2001 j l has shown a discrepancy between theoretical predictions 
and the observed main-sequence sample, especially at the region 
of the coolest stars. 

In this paper, we present a fiducial "single-star sequence" in 
the Hyades based on literature data and our own AstraLux Lucky 
Imaging observations. Further, we use the obtained sequence to 
provide a test for commonly-used stellar and atmospheric mod¬ 
els. The paper is organized as follows: Section 2 describes ob¬ 
servations and data from the literature; Section 3 discusses the 
theoretical models used; in Section 4 we show the "single-star 
sequence" in the Hyades compared to theoretical isochrones cal¬ 
culated using various stellar, and atmospheric models and deter¬ 
mine physical parameters of each Hyades "single-star"; in Sec¬ 
tion 5 we discuss our results and give a brief summary. 


tience et al. (|1998|l, Mermilliod et alT] ( |2009| l, |MorzinskI| ( |2011| l, 
Duchene et al. (2013|l and archival data from Hipparcos, Hub¬ 


ble Space Telescope (HST) and Washington Double Star (WDS; 
[Mason et al.|200l] ) catalogs. The HST data have also been pre¬ 
viously studied by Gizis & Reid ( [1995 1, Reid & Gizis \i991) 
and| Reid & Mahoney [([2000| l. Characteristics of the surveys are 
summarized in Table [iJ 


2.2. AstraLux lucky imaging observations 

The Lucky Imaging technique is based on series of short (few 
to several 10 ms) exposures, which "freeze" the speckle clouds 
caused by atmospheric turbulence. The series of frames is sorted 
by the best quality (based on the brightest pixel), and typically 
the best 1% to 10% of the images are shifted and co-added, re¬ 
sulting in a close to diffraction limited image. AstraLux Norte 
is a lucky imaging camera mounted on the 2.2m telescope at the 
Calar Alto observatory in Spain ( [Hormuth et al.|2008l l. The cam¬ 
era has a field of view of 24" x 24" (FOV) and a pixel scale of 
~ 47 mas/px. 

The observations were performed in Nov. 2011, Nov. and 
Dec. 2012. Depending on the brightness of the target and ob¬ 
serving conditions, either 15 or 30 ms exposure time was cho¬ 
sen, so that 20,000 or 10,000 short exposures were obtained for 
each target, respectively. Only one quadrant of the detector was 
read-out to facilitate shorter integration times. 

The Hyades targets for the AstraLux run are pre-selected 
based on the results of the literature and archive check. In to¬ 
tal 198 Hyades members are observed in SDSS i' and z' filters. 
The analysis for each target is done based on the best 1% frames 
that are co-added together. Possible companions are identified 
by a visual check of each co-added frame (see Fig. [T]) result¬ 
ing in 40 Hyades members that reveal one or several companion 
candidates. 28 of the 40 systems with companion candidates did 
not have any binary/multiplicity record in previously published 
surveys. Most of the AstraLux candidates to binary and multiple 
systems have only one epoch of observations, therefore they can¬ 
not be confirmed as common proper motion companions. Non¬ 
detections exclude companions outside a projected separation of 
3-7 AU (depending on a distance to a star), though these stars 
could still have unresolved companions closer-in. 


2. Observations and literature data 


To test interior and atmospheric models we use JHK^ 2MASS 
data ( Skrutskie et al.|2006 1 which allows us to get a homogenous 
photometric set for all the 724 possible members of the Hyades. 
The angular resolution of 2MASS is limited to ~ 3". To avoid 
a scatter on color-magnitude diagrams that can be introduced by 
unresolved binary and multiple systems, or optical blends with 
unrelated field or background stars, we check the literature and 
archived data bases for multiplicity and perform Lucky Imaging 
observations. The saturation threshold of the 2MASS photom¬ 
etry is ~ 4 mag, therefore fluxes for the brightest stars are 
replaced by the data from [Carney ( 1982[ ), with a transformation 
to the 2MASS system applied. The color uncertainties are cal¬ 
culated in the standard way exploiting independent photometric 
magnitude uncertainties estimates. 


2.1. Literature data 

We check previous studies that make an attempt to identify bi¬ 
nary and multiple systems in the Hyades including works ofH 
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In total 463 Hyades members are observed in one or more 
surveys and 213 of them are identified as candidates for binary or 
multiple systems leading to 250 stars that do not have identified 
companions outside projected separation of 2-7 AU (depending 
on survey and distance to the star). Importantly, the angular res¬ 
olution of instruments in the surveys we use varies from 0705 
to O'.'ll, which, additionally to various stars’ distances, gives 
a range of separation where companions can be detected (2-7 
AU; see Table 2). Hence, this means different detection thresh¬ 
olds for companions depending on in which survey(s) the system 
was observed. Since this complicates the statistical analysis for 
constraining the number of binary and multiple systems, a ro¬ 
bust estimate of the multiplicity properties of the Hyades will be 
a subject to a separate paper. 


3. Models 

In this section we describe the stellar evolutionary models and 
synthetic spectra used to produce theoretical isochrones for the 
Hyades and compare them with observed single-star sequence. 
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Survey 

Telescope 

Method 

Angular Resolution 

Observed 

Binary/Multiple 

Patience et al. (1998) 

5m Hale telescope 
Palomar observatory 

2.2 fim 

speckle imaging 

O'.'ll 

163 

33 

Mermilliod et al. (2009) 

Swiss Im telescope 
Haute-Provence Observatory 

CORAVEL 

spectroscopy 

- 

139 

25 

Morzinski (2011) 

Keck and Lick 
observatories 

Adaptive 

optics 

0'.'06 

75 

30 

Duchene et al. (2013) 

Keck II 

Adaptive 

optics 

0'.'06 

9 

5 

HST 

Hubble Space 
Telescope 

Imaging 

0'.'05 

57 

8 

Hipparcos catalogue 

Hipparcos Space 
Mission 

Imaging 

O'.'IO 

195 

22 

AstraLux Norte 

2.2m telescope 

Calar Alto observatory 

Lucky 

imaging 

O'.'ll 

198 

40 


Table 1; Characteristics of various multiplicity surveys in the Hyades. Some stars were studied by several surveys. The WDS catalog 
is not presented in the table because its instrument characteristics is inhomogeneous. 


HYA214 



1 arcsec 


Fig. 1: AstraLux Norte z'-band image of a binary candidate in 
the Hyades. The image is shown in log scale. 


3.1. Interior models 


We use two different stellar interior models - DARTMOUTH and 
new PISA grids produced for the Hyades case. Most recently, 
new BHAC15 ( Baraffe et al.|2015] l stellar models have been pub¬ 
lished. These models are adherent to the widely-used BCAH98 
models ( [Baraffe et al.|1998] l. However, BHAC15 models do not 


provide interior structures for super-solar metallicities yet, there¬ 
fore, we do not use them for the isochrone comparison. 

The Dartmouth Stellar Evolution Program (DSEP) ( [Potter 
[et al.||2008| ) is a family of stellar evolutionary models explor¬ 
ing a wide range of metallicities and a-enhancements. The basic 
equation of state for tracks M > O.SMq is a general ideal gas 
EOS with Debye-Hiickel correction ( [Chaboy^ & Kim||191^| l. 
The EreeEOS in the EOS4 configuration is used for the lower 
mass tracks. Evolutionary tracks are computed for masses from 
0.1 to 4 Mq. We use isochrones for [Ee/H] = 0.14 and [cr/H] = 
0.0. Eor the Hyades metallicity, the DSEP models include a con¬ 
vective core overshooting of 0.2 times pressure scale height for 
M> 1.1 Mq. 

The new PISA models have been computed by means of 
the ERANEC stellar evolutionary code (DeglTnnocenti et al. 
|2008|l, adopting input physics similar to those already discussed 


in |Tognelli et al.| ( [20li[ l and DeU’Omodarme et al.[ ( 2012] ). The 
main difference with respect to the models already available on 
the Pisa database pag^is the adoption of the SCVH95 EOS for 
the computation of stellar models with mass lower than 0.2 Mq. 
We adopt the recent [Asplund et al.[ ( [2009[ ) solar metal distribu¬ 
tion, and the corresponding mixing length parameter calibrated 
on the Sun, namely aML = 1-74. We also include a mild con¬ 
vective core overshooting (jSov = 0.2) for M > 1.2 Mq (Tognelli 
[et al.|201^ . The models have been computed in the mass range 
[0.1, 2.8] Mq from the early pre-MS evolution up to the ex¬ 
haustion of the central hydrogen. The corresponding isochrones 
in the age interval [400, 800] Myr, with an age spacing of 10 
Myr, have been generated. The stellar models have been com¬ 


puted for [Ee/H] = h- 0.14, which adopting the Asplund et al. 
(Z/2f)o = 0.0181 and AT/AZ = 2 (Casagrande [2007] l, corre- 
sponds to T = 0.283 and Z = 0.0175. 

To obtain a consistent set of magnitudes for both the Dart¬ 
mouth and Pisa models, we transform the theoretical isochrones 


'hrp://astro.df.unipi.it/stellar-models/ 
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from the (log Teff, logL/Lo) plane into the color-magnitude di¬ 
agram by means of our own calculation of photometric band. To 
do this we use synthetic spectra obtained from detailed atmo¬ 
spheric models following the method described in Section 


3.2. Atmospheric models 


To calculate synthetic photometry we use synthetic spectra pro¬ 
vided by BT-Settl atmospheric models ( [Allard et al.|2013] ). The 
BT-Settl models take into account gravitational dust settling in 
atmospheres of objects at temperatures below ~2600 K, follow¬ 
ing the approach described in |Rosso^ ( |1978| l. Opacities are in¬ 
troduced line-by-line to account for the effect of molecular ab¬ 
sorbers, as described in jAllard et al. (20031. The synthetic spec¬ 
tra are calculated using the radiative transfer model atmosphere 
code PHOENIX that implements static and radial (ID) approxi¬ 
mations ( Allard et al.|200l] l. We apply the latest publicly avail¬ 
able release of BT-Settl (last update in March 2015) that uses the 
Caffau et al. (20111 solar abundances. For effective temperatures 


Teff < 7000 K we use the BT-Settl 2010 version that exploits 
solar abundances ofjAsplund et al.|([2009|l. 


3.3. Synthetic photometry with TA-DA 


TA-DA stands for the Tool for Astrophysical Data Analysis (Da 
|Rio & Robberto|2012| l. TA-DA is an interactive software allow¬ 
ing to analyze stellar photometric data in comparison with theo¬ 
retical models and to derive stellar parameters using multi-band 
photometry. TA-DA is able to interpolate stellar interior models 
and to produce synthetic photometry by converting stellar pa¬ 
rameters into photometric magnitudes in given filters using grids 
of synthetic spectra. The conversion is done in the standard way, 
by integrating the synthetic spectra over the filter bandwidths 
and normalizing onto a spectrum of Vega. Therefore, TA-DA 
is well-suitable for combining various stellar interior and atmo¬ 
spheric models to derive photometric magnitudes for different 
evolutionary tracks and isochrones. We use stellar interior mod¬ 
els and synthetic spectra from grids of atmospheric models to 
obtain various isochrones for the Hyades with TA-DA. 


4. Hyades single-star sequence 

In this section, we give the Hyades single-star sequence as an ex¬ 
ample for testing theoretical models. Additionally, we determine 
stellar parameters of the sequence members by comparing pho¬ 
tometric observations to the calculated model isochrones. Abso¬ 
lute JHKs magnitudes for the stars on the single star sequences 
(see Table 2) are derived from the distance estimates in [Roserj 
|et al.| ( |2011| ) and the available near-infrared (NIR) photometry. 
For all stars we assume the NIR foreground extinction to be neg¬ 
ligible. 

4.1. Isochrone comparison 

We compute theoretical isochrones for combinations of BT-Settl 
2010 atmospheric models with PISA and DARTMOUTH inte¬ 
rior models for the age of 630 Myr, which is close to the 625+50 
Myr estimate of [Perryman et al.] ( |1998| l. We place the obtained 
isochrones on J vs. J - K^, J vs. J - H and Ks vs. H - color- 
magnitude diagrams together with the observed single-star se¬ 
quence of the Hyades. As can be seen, both PISA and DART¬ 
MOUTH tracks predict the observed sequence reasonably well, 
even the behavior of the "knee-shaped" part, around ~ Q.6Mq 
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which had proven to be problematic as previously has been 
shown by Roser et al.[ ( 201 l[ l for older generation evolution¬ 
ary models. The PISA models are available with (J3 - 0.2) and 
without convective core overshooting. Both sets of models are 
in good agreement with the presented data, the only difference 
being the cluster inferred age. As expected the inclusion of con¬ 
vective core overshooting leads to older ages. Adopting fi - 0.2 
the best fitting isochrone provides an age of 630 Myr, whereas 
isochrones without core overshooting give 550 Myr. 

The observed sequence reveal a larger scatter towards lower 
masses. This is because, due to the faintness of the objects, both 
photometric and kinematic measurements have larger errors. Ad¬ 
ditionally, a small (~0.05 mag) systematic offset can be noted on 
J vs. J - Ks and KsVS H - Ks diagrams - the observed sequence 
shifted to redder colors from the theoretical predictions. 


4.2. Physical parameters of single stars with TA-DA 


TA-DA (see Section [33] l includes an option that gives theoretical 
predictions for stellar parameters based on a comparison of syn¬ 
thetic photometry with observed photometric magnitudes. The 
dimension of the observational space (number of colors or mag¬ 
nitudes) must be greater or equal to the number of free parame¬ 
ters given by stellar interior models. The TA-DA parameter fil¬ 
ter is described in Da Rio & Robberto (2012 1 and performs a 
multi-band least-square fit. TA-DA also allows to estimate un¬ 
certainties for derived stellar parameters using a Monte Carlo 
simulation, in which the photometry is displaced according to 
photometric errors. 

Using TA-DA we derive theoretical stellar parameters for 
250 members of the Hyades single-star sequence (short elec¬ 
tronic version of table Table 2) using BT-Settl2010+PISA 
isochrones for [Fe/H]=+0.14. 


4.3. Spectral Energy distribution from X-ray to MIR 

For all stars of the single star sequence, we compiled the avail¬ 
able literature data from the X-rays to mid-infrared ROSAT X- 
ray ([Stern et al.[[T995l), GAFEX UV ([Martin & GAFEX Team| 


|2005| , APASS gri ( [Henden et al.|[20l5[ ), 2MASS JHKs, and 

WISE W1 to W4), and computed absolute fluxes and magnitudes 
based on the parallaxes from Roser et al. ( 2011[ l. The values are 
provided in the long electronic version of Table 2. 


5. Summary and future prospects 

Using previously published high angular resolution and spectro¬ 
scopic studies and our own Fucky Imaging AstraFux Norte ob¬ 
servations, we compiled a "single-star sequence" for the Hyades, 
which at an average distance of 45 pc is the open cluster closest 
to the Sun. In total, we identified 250 members of the cluster 
that do not show signatures of a companion outside 2-7 AU pro¬ 
jected separation (depending on the distance to each individual 
member and the survey). 

Comparison of the near-infrared properties of the single-star 
sequence with theoretical isochrones based on PISA and DART¬ 
MOUTH stellar interior models and BT-Settl 2010 atmospheric 
models shows an overall good agreement for the mass range 0.13 
to 2.6 Mo. The only disagreement between models and observa¬ 
tions is that for masses below 0.6 Msun the observed J - Ks and 
H - Ks colors are systematically redder than the isochrones by 
0.05 mag. But despite this, the isochrones calculated using most 
recent interior and atmospheric models, show a significantly bet- 
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ter agreement with observations than isochrones produced by 
older generation models. The improvement can be particularly 
seen around 0.6-0.8 Mo, where new isochrones are able to re¬ 
produce the "knee-shaped" part of the observed sequence. 

We also estimate physical parameters (mass,effective tem¬ 
perature and surface gravity) for the Hyades single-star sequence 
members based on theoretical predictions of combined PISA in¬ 
terior and BT-Settl 2010 atmospheric models. 

The resulting "single-star sequence" can be used for test¬ 
ing various theoretical models and also for selecting candidates 
to search for close, previously unresolved binaries using spec¬ 
troscopy and high-resolution imaging instruments. In the near 
future, results from the Gaia survey will provide us with kine¬ 
matic, photometric and spectroscopic information that will help 
to reduce the scatter at the lowest-mass part of the observed 
sequence, which will enable us to test theoretical models with 
stronger confidence. 
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Fig. 2: J vs J - Ks (top left panel), J vs. J - H (top right panel) 
and Ks vs. H - (bottom right panel) color-magnitude dia¬ 
grams. The Hyades observed single-star sequence is over-plotted 
with theoretical isochrones (630 Myr) from different evolution¬ 
ary models - PISA (dashed blue line) and DARTMOUTH (solid 
red line). Synthetic photometry for isochrones are calculated in 
TA-DA using BT-Settl 2010 synthetic spectra. 
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Table 2: Absolute JHKs(2MASS) photometry and physical parameters of the Hyades single-star sequence members. The Resolution column correspond to values outside which 
no companion was detected. Objects observed only by the spectroscopic survey of Mermilliod et al. (2009 have no Resolution value. Hyades IDs are from Roser et al. ([201 1 . 
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